







1. Report No. 2. Government Accession No. 

NASA TM X-2636 

3. Recipient's Catalog No. 

4. Title and Subtitle 

SIMPLIFIED THERMOCHEMISTRY OF OXYGEN IN LITHIUM 
AND SODIUM FOR LIQUID METAL COOLING SYSTEMS 

5. Report Date 

> September 1972 

6. Performing Organization Code 

7. Author(s) 

Leonard K. Tower 

8. Performing Organization Report No. 

E-7021 

10. Work Unit No. 

503-25 

9. Performing Organization Name and Address 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland,. Ohio 44135 

11. Contract or Grant No. 

13. Type of Report and Period Covered 
Technical Memorandum 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D. C. 20546 

14. Sponsoring Agency Code 


15. Supplementary Notes 


16. Abstract 

Plots of oxygen chemical potential against composition of lithium - oxygen solutions and sodium - 
oxygen solutions for a range of temperature were constructed. For each liquid metal two such 
plots were prepared. For one plot ideal solution behavior was assumed. For the other plot, 
existing solubility limit data for oxygen in the liquid metal were used to determine a first-order 
term for departure from ideality. The use of the plots in evaluating the oxygen gettering capa- 
bility of refractory metals in liquid metal cooling systems is illustrated by a simple example 


involving lithium, oxygen, and hafnium. 


18. Distribution Statement 

Unclassified - unlimited 


19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21. No. of Pages 

22. Price* 

Unclassified 

Unclassified 

22 

$3.00 


17. Key Words (Suggested by Author(s) ) 

Lithium - oxygen solution 
Sodium - oxygen solution 
Liquid metal coolants 


For sale by the National Technical Information Service, Springfield, Virginia 22151 






















SIMPLIFIED THERMOCHEMISTRY OF OXYGEN IN LITHIUM 
AND SODIUM FOR LIQUID METAL COOLING SYSTEMS 
by Leonard K. Tower 
Lewis Research Center 

SUMMARY 

Plots of oxygen chemical potential against composition of lithium - oxygen solutions 
and sodium - oxygen solutions for a range of temperature were constructed. For each liq- 
uid metal two such plots were prepared. For one plot ideal solution behavior was as- 
sumed. For the other plot, existing solubility limit data for oxygen in the liquid metal 
were used to determine a first-order term for departure from ideality. 

The use of the plots in evaluating the oxygen gettering capability of refractory metals 
in liquid metal cooling systems is illustrated by a simple example involving lithium, ox- 
ygen, and hafnium. 


INTRODUCTION 

The use of liquid alkali metals as high temperature heat transfer media is being ex- 
tensively investigated. Familiar applications include forced convection cooling loops in 
nuclear reactors. Recently, heat pipes containing alkali metals have been proposed as a 
means of removing reactor heat (ref. 1). 

Because of the elevated temperature in these systems and the high chemical reac- 
tivity, some rather difficult materials problems are encountered. For instance, corro- 
sion of containment vessels by the liquid is observed. Solution of the vessel wall by the 
liquid metal has been studied (ref. 2) and found to be aggravated by small amounts 
of oxygen dissolved in the wall material (ref. 3). The formation of a soluble ternary 
oxide of the refractory and alkali metals is suspected as part of the corrosion mechanism 
(refs. 4 and 5) for some refractory metals. 

Because oxygen in liquid metals at elevated temperatures is apparently a cause of 
corrosion, various schemes have been proposed for alleviating the problem. The alloy- 
ing of oxygen -gettering elements with the metals used in containment vessel walls 
(ref. 6) and the use of gettering traps have been proposed. 



The success of such methods depends on the relative chemical potentials of oxygen 
in the liquid metal and in the gettering materials. An estimation of the chemical poten- 
tial of oxygen in the liquid metal becomes necessary in assessing the feasibility of getter- 
ing methods. 

No direct measurement of the chemical potentials of oxygen in the liquid metals is 
known to exist. This information must therefore be inferred from solubility limit data 
for oxygen in the liquid metals. This task is undertaken in the present report for lithium 
and sodium, two coolants of considerable interest for nuclear reactor application. 

The present report contains plots of chemical potentials for oxygen dissolved in lith- 
ium and in sodium as functions of temperature and concentration of oxygen in solution. 

For each liquid metal two such plots were prepared. The assumption of ideal solution 
behavior was made in constructing one plot. For the other plot existing oxygen solubility 
limit data was used to determine the coefficient of a first order term for departure from 
ideality. The form of the expression used for this nonideal correction was essentially 
that of the regular solution (ref. 7). The use of the chemical potential diagrams in as- 
sessing chemical stability of a system containing refractory metal, alkali metal and ox- 
ygen is illustrated by an example involving hafnium. 

Because of the uncertainties in the solubility limit data analyzed and the arbitrary 
nature of some of the assumptions, thermochemical data herein for the nonideal case is 
at best very tentative. The ideal case has been included for comparison so that the sen- 
sitivity of the computed thermochemical results to the assumptions made in the analysis 
can be assessed. The potential importance of liquid metal cooling systems to nuclear 
reactor technology makes such tools as devised here for the analysis of thermochemical 
problems necessary until more definitive information is available. 


ANALYSIS AND DISCUS.SION 

Coordinates of temperature - composition phase diagrams have sometimes been used 
to determine the thermochemical properties of solution phases. Various analytical meth- 
ods have been devised for this task, some of them suitable for binary metallic systems 
(refs. 8 to 10) and some for multicomponent systems (refs. 11 and 12). The methods 
consist of determining the values of constants in Taylor series or similar expansions re- 
lating the chemical potentials of the solution components to composite -variables. 

In a study (ref. 13) related to the present one, the thermochemical properties of a 
system containing cesium and oxygen was computed from an experimental phase diagram 
(ref. 14). This diagram covered a far greater composition range than the solubility limit 
data available for oxygen in lithium and in sodium. The method of references 8 to 10 as 
used in reference 13 is not suitable for oxygen in lithium and in sodium. 
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A very elementary approach will be employed herein. The first assumption made is 
that the components of the lithium - oxygen and sodium - oxygen systems are Li, LigO 
and Na, NagO. Experimental evidence shows the presence of O'" ions in liquid alkali 
metal systems containing oxygen (refs. 15 to 17). Even a small concentration of such 
ions might affect the electrical conductivity (ref. 17) or play a role in the kinetics of the 
corrosion mechanism (refs. 15 and 18). The absence of any quantitative data for concen- 
trations of O' " relative to total oxygen in liquid lithium or sodium precludes considering 
the effect here. 

The second simplifying assumption is the applicability of the so-called regular solu- 
tion approximation in a modified form. Regular solution behavior has been found in 
many binary systems involving metallic species (ref. 19). In systems exhibiting this be- 
havior, ju a and of species a and b are given by 

11^ = lil + RT \n(l -x) + Lx 2 ( 1 ) 

p b = + RT ln x + L (* ~ x)2 (2) 

where p° designates the free energy of the pure liquid at the temperature T in ques- 
tion, x is the atom fraction of b in solution, and R is the gas constant. Ideal solution 
behavior occurs when the interaction parameter L is zero. 

Sometimes nonideality in solutions is handled by defining activity coefficients y„ 

d 

and y^ such that (ref. 20) 


M a = + RT ln y a y 


(3) 


M b = MC + RT ln V 


(4) 


where x = 1 - y. The parameters y a and y^ generally vary with x. If y a — constant 
as y ~ 0, Henry’s law is said to be obeyed in the limit. If y — 1 as y — 1, Raoult’s law 
for ideal behavior is obeyed in the limit. Similar considerations apply to y^. The be- 
havior of equations (1) and (2) in the foregoing limits is seen to be nearly the same as 
that of equations (3) and (4). The principle difference is that L either is assumed con- 
stant or has a temperature dependence of the form L = L Q + L^T (ref. 19). 

The modification adopted here consists of setting L Q = 0. Equations (1) and (2) then 
become 


M a = M a + RT find - *) + Bx 2 ] 


(5) 
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( 6 ) 


M b = /xg + RT [in x + B(1 - x) 2 ] 

This modification will be shown to be advantageous in handling the experimental data for 
oxygen in lithium. 


Oxygen in Lithium 


Measurements of the solubility of oxygen in lithium are reported in reference 21. 

Solid dilithium oxide was placed in contact with lithium. The amount of oxygen remaining 
in the lithium after filtration was then determined by analysis. Table I presents the re- 
sults from reference 21. 

The oxygen concentrations are also converted in table I to equivalent mole fractions 
x of dilithium oxide (LigO) in a solution of LigO and lithium. This is consistent with the 
assumption herein that all dissolved oxygen is in LigO. 

At the condition of saturation, LigO in the solid phase must be in equilibrium with 
LigO in solution, and the free energy 0(c) of crystalline LigO must equal the 

chemical potential of dissolved LigO. 2 

Equations (2) and (6) for b as LigO were evaluated for parameters L and B from 
the x data in table I. Free energies of fusion, q^ c ^ , for LigO from refer - 

ence 22 used in the computations represent an extrapolation into the supercooled region 
for LigOQ.) since the temperatures are below the freezing point. The resulting L from 
equation (2) and B from equation (6) are listed in table I. 

The mean value of B for the data in table I is 1. 835 with a standard deviation of 
0. 221. No single satisfactory mean value of L can be determined from the data. If L 
is permitted to vary linearly with temperature, a much better fit can be achieved, the 
correlation being L = -1950 + 6.96 T. At the melting point of LigO, 1843 K (refs. 17 
and 19), the value of L is about 9970. 

Allowable values of L and B are limited by theoretical considerations discussed 
in the appendix. The maximum allowable value of B is shown there to be 2. The max- 
imum allowable value of L at the melting point is about 7300. Since the value of L ex- 
trapolated from the experimental data exceeds this, the B formulation of equations (5) 
and (6) was used herein with B = 1. 835. 

The solubility limit x gat for LigO in lithium computed from equation (6) using this 

B with (AF°) Li from reference 22 is shown in figure 1 against reciprocal temper - 

2 

ature. Also shown are the experimental data points from reference 21 used in determin- 
ing B. The ideal solubility limit determined from equation (6) by setting B = 0 is shown. 
In the temperature range of the experimental data about two-thirds of a decade separates 
the ideal and nonideal solubilities. 
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Equations (5) and (6) with B = 1. 835 together describe the thermochemical prop- 
erties of the Li - LigO system. For convenience in studies involving the presence of ox- 
idizable materials, the thermochemical base is changed from Li 2 0(s) and Li(s) to an at- 
mosphere of diatomic oxygen in equilibrium with the Li - LigO solution. For the reac- 
tion 


Li 2 0(s) Z 2Li(s) + I O a (g) 
2 


thermochemical equilibrium requires that 


(7) 


^Li 2 0(s) " 2fi Li(s) 


+ \ 2 & 


( 8 ) 


Here, 


M o„(g) ~ ( f t)o 0 ~ ( f t)o 9 + RT ln p O r 


(9) 


where the pressure p for thermochemical purposes is in atmospheres (1 atm = 
1.013x10 s N/m). Combination of equations (5), (6), and (9) yields 


ln 


^AFn 


^AF, 


p O 




2 \RT/0 2 \RT/Li 2 0(c) \RT /LigOfc,!) 


+ ln 


+ B(1 - 2x - xj (10) 


(1 - xY 


where the standard state free energy of formation of Li 2 0(c) is 


( AF f )Li 9 Q(c) " ( F T)Li 9 0(c) ' 2 ( F t)l1(1) - 9 ( F t)o 


(ID 


Values of [(F^,) - (H 29 g)]/ T l° r LLjO(c) and Li(c,l) were taken from reference 22, while 
K F t) - H 298)] /T for °2 (g) and ( H f,298)Li f 0(c) were l rom reference 23. A plot of equa- 
tion (10) appears as figure 2(a) herein. Figure 2(a) shows (aF^/Rt)q against mole 
fraction of LigO in solution, for several temperatures. 

By setting B = 0 equation (10) was also evaluated for the ideal case. The results 
are shown as figure 2(b). Comparison of figures 2(a) and (b) reveals that (AF° /RT)q 

2 

is higher (less negative) by about 4 units in the nonideal case. This amounts to about 
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23 000 joules per mole (5500 cal/mole) at 700 K and about 59 000 joules per mole (14 000 
cal/mole) at 1800 K. 


Oxygen in Sodium 

Many determinations have been made of the solubility of oxygen in sodium. Refer- • 
ence 16 summarizes and discusses data of many investigators taken prior to 1966. This 
body of data together with subsequent information was evaluated critically in refer- 
ence 24. The following standard equation for the solubility of oxygen in sodium was rec- 
ommended as a result of the evaluation: 


log 1Q s = 6. 239 - 


2447 


(12) 


where s is weight parts per million of oxygen in sodium. Converted to mole fraction x 
of sodium oxide in solution equation (12) gives 


log x = 14. 37 - 


5634 


(13) 


By using equation (13) together with equation (6) and values of (AF°) Na q^ c ^ from ref- 

erence 22, B was evaluated at several temperatures encompassing the range covered by 
the many investigators whose data contributed to the construction of equation (12). These 
temperatures and the corresponding values of B are 400 K, 1.75; 500 K, 1.92; 600 K, 
2.13; 700 K, 2.33; 800 K, 2.53. The mean of these values of B is 2. 13. 

While B is understood to be a constant, the values determined from the curve fit of 
reference 24 show a decided trend with temperature. Another approach would have been 
to solve for B for each data point used in constructing equation (12) and to take the 
mean. The amount of scatter in the data together with the number of assumptions in- 
volved in the present analysis makes such an effort unwarranted. 

The mean B of 2. 13 deduced from equation (12) exceeds the maximum allowable 
value of 2 discussed in the appendix. For the Na - Na 2 0 studies herein a value of B = 2 
will be adopted. 

Figure 3 shows the solubility limit x gat against reciprocal temperature for NagO in 
sodium, computed from equation (6) using a B of 2 and (AF°) Na q^, ^ from refer - 

ence 22. The ideal solubility limit found by setting B = 0 in equation (6) is also shown, 
together with equation (13) from reference 16. 
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Figure 4(a) is a plot of (aF^/RT)q against mole fraction of NagO in solution with 

sodium at several temperatures for oxygen in equilibrium with the solution. The data 
were computed for B = 2, the nonideal case, using equation (10), with thermochemical 
data from reference 22. Figure 4(b) is the corresponding plot for B = 0, the ideal case. 


Example Using Oxygen-Contaminated Lithium in Contact with Hafnium 

The use of figures 2 and 4 in evaluating the chemical stability of a system containing 
a refractory metal in contact with oxygen-contaminated lithium can be illustrated by an 
example involving hafnium. Hafnium is a common constituent of alloys that are being 
considered for structures containing liquid alkali metal coolants. The hafnium in part is 
intended as a getter to counter oxidation of the other alloying constituents. 

The equilibrium between Hf0 2 (c) and Li 2 0(s) is expressed as 

2Li 2 0(s) + Hf(c) =Hf0 2 (c) + 4Li(s) (14) 


which requires that 


2 ^Li 2 0(s) + ( F T)Hf(c) ~ ( F T)HfQ 2 (c) + 4 ^Li(s) 


(15) 


Equilibrium (14) can be expressed as the sum of the following equilibria: 


2Li 2 0(s) 3t4Li(s) + Og(g) 

(16) 

Hf(c) + 0 2 (g) rHf0 2 (c) 

(17) 


The condition for equilibrium (16) to exist is given by equation (8), which can be reexpres- 
sed as 


2fx Li 2 0(s) ' 4 ^Li(s) " ( F T)o 2 (g) “ (AF T ) 0 2 (g) 
For equilibrium (17) to exist, 

( F T)Hf(c) + ^T^g) = ( F T)Hf0 2 (c) 


(18) 


(19) 
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Subtracting 


( F T)o 2 (g) 


from both sides of equation (19) gives 


( AF T)Hf0 2 (c) " ( AF T)o 2 (g) 


( 20 ) 


where (^ F f)jifo *- s the standard state free energy of formation for Hf0 2 (c). When 
2 

(aFj^q ^ in equations (18) and (20) are equal, these equations can be combined to ob- 


(g) 

behavior , 

and on figure 2(b) for ideal solution behavior. The HfOg data are from reference 25. 

The line and its x coordinate will henceforth be designated as x(Hf0 2 ). 

For NagO in solution (fig. 4) the line for x(Hf0 2 ) cannot be plotted since its small 
magnitude places it too far to the left. The effectiveness of hafnium in gettering oxygen 
from sodium is thus illustrated to be much greater than in gettering oxygen from lithium. 

For LigO concentrations to the left of the line x(Hf0 2 ) on figure 2, Hf0 2 (c) will not be 
stable, while to the right of the line it will be. To determine on which side of the line 
x(HfC> 2 ) a given Li - Hf - O system lies, an inventory must be made of the total moles of 
lithium (N Li ), of hafnium (N Hf ), and of atomic oxygen (Nq) in the system. All phases 
and all compounds must be considered in the inventory. For this example the gas volume 
will be assumed small enough to contain negligible material. 

A simple criterion for the stable existence of Hf0 2 (c) can be readily derived. In the 
presence of stable Hf0 2 (c) the excess oxygen forms LigO, the number of moles of which 
are 


tain equation (15), the condition for equilibrium (14). A line of (AF t /RT)q 

2 

(AF°/RT) Hf Q ^ (see eq. (20)) is plotted on figure 2(a) for nonideal solution 


N Li 2 G = N 0 - 2N Hf 


(21) 


Of the total amount of lithium, in the system, the amount as free lithium is 


N Li ~ N Li “ 2N Li 2 0 


( 22 ) 


The number of moles of solution is 


N Sol = N ’li + N Li 2 Q 


( 23 ) 


The mole fraction x of LigO in solution is then 




8 



(24) 


x = 


"Li 2 ° 


N, 


Sol 


N 0 ~ 2N Hf 
N Li + 2N Hf " N 0 


For Hf0 2 (c) to be thermochemically stable 

x > x(Hf0 2 ) 


or 


N 0 - 2N Hf 
N Li + 2N Hf - N 0 


> x(HfQ 2 ) 


(25) 


If criterion (25) is not satisfied, Hf0 2 (c) cannot be stable. In such a case the concentra- 
tion x will be less than x(Hf0 2 ). 

The example for HfOg can be extended to other hafnium oxides and to solutions of ox- 
ygen in hafnium when currently unavailable thermochemical data become available. Re- 
cently the compounds Hf^CXc) and HfgO(c) have been identified (ref. 26). Lines for these 
compounds should lie to the left of Hf0 2 (c) in figure 2. Solid hafnium can contain up to 
7 percent oxygen in solution (ref. 26). The lines of (AF T /RT)C> 2 for varying concentra- 
tions of dissolved oxygen in hafnium will lie still further to the left of x(HfC> 2 ) in figure 2. 

Hafnium as a minority element in some refractory metal alloys such as T- 111 
(tantalum - 8 tungsten - 2 hafnium) is in a single phase solution. For this case ( F ^)Hf(c) 
in equation (15) must be replaced by a chemical potential pertaining to hafnium in the 
metal solution. An additional complication is the possible formation and solution of 
LigHfOg in lithium (ref. 2). 

The foregoing example, despite the dearth of thermochemical data for hafnium ox- 
ides and for oxygenated lithium, illustrates the procedures involved in estimating the ef- 
fectiveness of getters for removing oxygen from alkali metals. The constructions in fig- 
ure 2 provide a first order estimate of this gettering effectiveness for the Hf - O - Li 
system. 


CONCLUDING REMARKS 

Plots of oxygen chemical potential against composition were presented for solutions 
of oxygen in lithium and in sodium at several temperatures. Because of the uncertainties 
in the experimental data used to construct the diagrams for nonideal solution behavior 
and the assumptions required, plots were also presented for ideal solution behavior. 
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Comparison of the results for the ideal and nonideal cases gives an indication of the 
sensitivity of the thermochemical results to the assumptions made in the nonideal anal- 
ysis. The potential importance of liquid alkali metals in nuclear reactor cooling systems 
justifies the use of the tentative data presented herein in analyzing problems involving 
oxygen contamination until better data can be derived. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 7, 1972, 

502-25. 
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APPENDIX - UPPER L1AA1TS OF L AND B 


Upper limits of L and B are imposed by thermochemical requirements. At fixed 
temperature 

~ ( F T)b(c) ^ Al) 


for x < x sa f Otherwise, the purported x ga ^. must be conceded to be erroneous and the 
saturation line at that temperature shifted to a smaller value of x where does in- 
deed equal ( F T)b(c)‘ Fi S ure 5(a) illustrates the situation which would exist in such a 
case. This figure is drawn in such a manner as to suggest that one requirement for the 
saturation point x a j. is that the variation of ju^ with x be monotonic at x < x„ 

For equation (2) to be monotonic, 


sat‘ 



— - 2(1 - x)L > 0 
x 


(A 2) 


for xSx sar 

A possible variation in which satisfies inequality (Al) but not inequality (A2) is 
sketched in figure 5(b). The situation depicted requires that a point of inflection lie be- 
tween two real roots of inequality (A2) . For a point of inflection 


= . 31 + 2L = 0 

\Sx 2 /T x 2 


the only positive abscissa of which is 


x = 



(A4) 


The locations of the maxima or minima of inequality (A2), where (3ju b /9x) T = 0, are 

x=ia JT. 

2 2 f 


2RT 


(A 5) 


If the point of inflection is to lie between the roots of 0p^/3x as depicted, equations (A4) 
and (A5) require that 
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and 



(A9) 


Squaring again and rearranging result in 


il > 2 (A10) 

RT 

However, for real x, equation (A5) requires that L/RT < 2. Except for the special case 
L/RT = 2, x = 1/2 these are contradictory requirements. The conclusion follows that the 
situation in figure 5(b) cannot exist for chemical potentials obeying equation (2). The var- 
iation of with x must therefore be monotonic, with inequality (A2) applying for 

x “ x sat' 

From inequality (A2) 


T.< ^ 

2x(l - x) 


(All) 


The minimum value of the right side, occurring at x = 1/2, is 2RT. Therefore, the 
maximum allowed value of L is 2RT. 

Precisely the same arguments, applied to equation (6), give B < 2. 
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TABLE I. - ANALYSIS OF SOLUBILITY DATA FOR 


OXYGEN IN LITHIUM FROM REFERENCE 21 


Temperature, 

°C 

Oxygen in 
lithium , 
wppm 

Mole fraction of 
LigO in solution, 
X 

L 

B 

250 

58 

2 . 516 X 10' 6 

2185 

2 . 102 


74 

3.210 

1932 

1.858 


85 

3.687 

1788 

1.720 


92 

3.990 

1706 

1.641 


96 

4 . 164 

1662 

1.598 


101 

4.381 

1609 
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Figure 2. - Chemical potential diagram for diatomic oxygen in equilibrium with a solution of lithium and l_i 2 0. 
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Figure 2. - Concluded. 
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(a! Nonideal case, B = 2.0. 

Figure 4. - Chemical potential diagram for diatomic oxygen in equilibrium with a solution of sodium and Na 2 0. 
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